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Due to increasing research being conducted at ICASE in the field of fluid mechanics, 
future ICASE reports in this area of research will be printed with a green cover. Applied 
and numerical mathematics reports will have the familiar blue cover, while computer science 
reports will have yellow covers. In all other aspects the reports will remain the same; in 
particular, they will continue to be submitted to the appropriate journals or conferences for 
formal publication. 
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THE THREE-DIMENSIONAL FLOW PAST A RAPIDLY 
ROTATING CIRCULAR CYLINDER 


James P. Denier 1 
School of Mathematics 
University of New South Wales 
AUSTRALIA 
and 

Peter W. Duck 1 
Department of Mathematics 
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UNITED KINGDOM 


ABSTRACT 

The high Reynolds number (Re) flow past a rapidly rotating circular cylinder is inves- 
tigated. The rotation rate of the cylinder is allowed to vary (slightly) along the axis of the 
cylinder, thereby provoking three-dimensional flow disturbances, which are shown to involve 
relatively massive (O(Re)) velocity perturbations to the flow away from the cylinder surface. 

Additionally, three integral conditions, analogous to the single condition determined in 
two dimensions by Batchelor (1), are derived, based on the condition of periodicity in the 
azimuthal direction. 
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supported by the Australian Research Council and by SERC. 
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1 Introduction 


One of the most important and fundamental theorems of Fluid Mechanics to have been 
developed in relatively recent times is the closed streamline theory of Batchelor (1). This 
states that in a two-dimensional flow involving closed streamlines 

<j> curl u> ■ ds = 0, (1) 

where u> is the vorticity and ds the line element along a streamline. At large Reynolds 
numbers, when viscous effects can be neglected, this condition, taken with the inviscid 
(Euler) equations of motion, leads to the result that the vorticity is constant within a region 
of closed streamlines. 

This theorem was generalised by Blennerhassett (2) to the situation which, rather than 
being one of closed streamlines, involves helical streamlines. Specifically a three-dimensional 
flow with three velocity components, all of which are independent of the axial direction; this 
implies constant axial pressure gradient. This analysis, in addition to confirming Batchelor’s 
(1) result, leads to an additional result linking the constant axial pressure gradient and the 
viscous terms. This result, taken in conjunction with the inviscid flow equations yields the 
result that the axial velocity must be proportional to the streamfunction for the motion in 
the plane normal to the axial vector. 

Little progress appears to have been made in further extending the ideas of Batchelor 
(1) into the three-dimensional regime, in spite of the importance of flows of this type. It is 
the aim of this paper to address this issue. 

Two-dimensional flow's involving closed streamlines involve additional interesting sub- 
tleties. Even though the work of Batchelor (1) involves the inclusion of (small) viscous 
effects, there is insufficient flow physics in the constant vorticity result to actually determine 
its value. This must be determined by recourse to boundary-layer regions located on body 
surfaces. Riley (3) presented an example of two-dimensional flow inside an elliptic container, 
driven by slippage of the container w r alls. By Batchelor’s (1) theorem the flow in the core of 
the container must be that of uniform vorticity, and Riley (3) showed how just one particular 
value of the vorticity produced the appropriate behaviour (exponential decay at the outer 
edge) inside the wall boundary layers. These ideas have subsequently been extended by other 
authors to other situations (see below), including situations involving tw r o w r all layers, see 
Duck (4). 

One of the classical results of potential-flow theory is that of uniform flow past a circular 
cylinder, with superimposed circulation. Physically this circulation may be caused by the 
rotation of the cylinder. A number of authors have studied the relationship between the rate 
of rotation and this circulation. Glauert (5) considered the large Reynolds number, large 
rotation rate problem, whilst Moore (6) considered the finite Reynolds number, large rotation 
rate problem. Loc. (7) and Ingham (8) have considered fully numerical solutions (finite 
Reynolds numbers, finite rotation rates), whilst Nikolayev (9), Negoda and Sychev (10), 
Sychev (11), Lam (12) have investigated the large Reynolds number (finite rotation rate) 
problem, partly using the ideas of Riley (3). There appears to be a monotonic relationship 
between rotation rate and circulation, and the simple model of outer, potential flow together 
with a unidirectional boundary layer is only appropriate for rotation rates above some critical 
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value. As this critical value is approached, a stagnation point forms off the surface of the 
cylinder wall, but inside the boundary layer; this is sufficient to disrupt the entire model at 
lower values of rotation rate. 

In this paper we consider a three-dimensional analogue of the above problem. We take 
a uniform, straight circular cylinder, and a uniform flow far from the cylinder directed 
perpendicular to the axis of the cylinder, and the cylinder surface is rotating with a large 
angular velocity, which is dependent on axial location. We further suppose that the Reynolds 
number (as defined in Section 2) is large. The restriction on large angular velocity renders 
the problem tractable to analytic techniques, whilst pointing the way in which more general 
rotation rates may be tackled. This analysis is presented in Section 3, and the corresponding 
numerical results are in Section 4. In Section 5 we formulate three integral conditions, 
analogous to the condition of Batchelor (1) described above. In Section 6 we present our 
conclusions. 
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2 Formulation 


We take a straight circular cylinder of radius a, together with polar coordinates ( ai\9,az ), 
with r measured radially, 9 azimuthally in an anti-clockwise direction and 2 axially. The 
surface of the cylinder is rotating with angular velocity uj + 7 cos Aar, where u, 7 and A are 
constants. Referred to the coordinate system described above, we suppose that far from the 
cylinder the fluid velocity takes the form (/^(cos 9, — sin 8, 0). The flow Reynolds number is 
defined by 

Re = (2) 

V 

where v is the kinematic viscosity of the fluid (assumed constant). In this paper we are 
primarily interested in the regime Re 1. 

Two further non-dimensional parameters may be defined, namely 


n = 


uja 

uZ' 


( 3 ) 


and 



( 4 ) 


describing the rate of rotation and degree of three-dimensionality introduced into the prob- 
lem. 

To make further progress in this paper, we assume that |e| <C 1 and further that |fl| 1 
(but that e = o(fi -1 )). Thus we focus our attention on small amounts of three-dimensionality 
and high rotation rates. 

We write the velocity vector as UooU, and the pressure as pU^p, where p is the density 
of the fluid (assumed constant). We then develop the solution in powers of e, namely 


u = u 0 + cu + 0(e 2 ), (5) 

p = Po + cp + 0(e 2 ). ( 6 ) 

For r — 1 = 0(1), when Re > 1 we assume the basic flow takes the form 

u 0 = (cc^l-i), — sin«(l + i) + 5o), (7) 

this being merely the potential flow solution. There also exists a boundary layer wherein 

?• — 1 = 0(Re _1,/2 ). The problem for Uo in the boundary layer has been studied in the past 

by a number of authors, as noted in the previous section, however as ft — + 00, 


u 0 — > (0,0,0) + 0(1) 


( 8 ) 


for r — 1 = 0(Re 1 ) . 

The full problem for u is 


V • u = 0, 


( 9 ) 


with the following boundary conditions 


( 11 ) 

(12) 


u = (0,cosAz,0) on r = 1, 

u — ► 0 as r — > oo. 

In the following section we consider the solution of (10) as and Re — » oo. 
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3 The limits Re — » oo, Q — > oo 


Since the perturbation velocity u is triggered by an 0( 1) amount (see (11)) it is tempting 
to speculate that this in turn will lead to general perturbation quantities away from the 
cylinder of 0(1). However this is certainly not the case. It turns out that the only plausible 
solution for r = 0(1), i.e. the only asymptotic development that leads to self-consistency 
involves magnitudes very much larger than 0(1). 

Considering first the solution for r — 1 = 0(1), then on account of the linearity of the 
system (9), (10), we may write 

u = (u(r, 0) cos Az, u(?’, 0) cos Az, u>(r, 0) sin Az), (13) 

p = p(r, 0)cosAz. (14) 

We also have (see (7)) that as Re -» oo, 

uq — ► (Voi,fif/oo + f7 O i,0), (15) 


where 



U m = — (l + sin 0. (16) 

The only meaningful series development of the perturbation velocities turns out to be 


v = fl Re u 0 (r) + ■ • • + Re tq(r, 0) 4- ■ • • , 

u = Re uj (r, 0) -i f- u 2 (r) 4 , 

w = 0 Re w 0 (r) 4- • • • + Reiei(r,0) + • • • , 

p = ft Rep o (*',0) + • ' 1 + ftpi(r) H ■ (17) 


These expansions prove to be the key to the solution of our problem, although a priori , it 
is difficult to justify these expansions, and so a posteriori verification is therefore necessary. 
In Section 5, some prima facie justification for this form of solution development is given. 

Consider now the leading order terms in (17) that are independent of 0 . It turns out, 
perhaps surprisingly, that the flow is not governed by the inviscid Euler equations, but 
rather is predominantly viscous in nature. Taking O(fl) terms in the radial momentum 
equation yields the following equation 


2 u-2 dp x 1 2 v Q 

— J b Vorr 4 VOr ~ O V 0 

r i r r i 


The 0(Re } ) terms in the azimuthal momentum equation yield 


. 1 \ 2 u 'i _ n 

U 2 rr 4“ ^‘2r ^ ^'2 .5 6) 


(18) 


(19) 
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whilst the 0(17) terms in the axial momentum equation yield 

U^Orr 4” Or A WJ() — Aj^| , (20) 

The continuity equation, to leading order, may be written 

Vor + -Vo + A Wo = 0. (21) 

r 

With this system there is no difficulty in imposing the full and correct boundary condi- 
tions; we therefore have 


M r=1 ) = w o(r = 1) = 0, 

u 2 (r=l) = 1, (22) 


with 

no, Wo, \i 2 — » 0 as r — » oo. 

The solution for u 2 can be written in terms of Bessel functions, namely 

_ ff,(Ar) 

“ 2 *i(A) ’ 

whilst it is possible to eliminate w o and ;q, from (18), (20), (21) to yield 

2 hi ( 3 


tut 2 m ( O % tf 
v o + ~ v o - + 2AM v 0 


3 2A 2 


+ I “T ~ 


2A 2 


,3 


7" r° + i A = 


2tx,A 2 


(23) 


(24) 


(25) 


to which the following additional (alternative) boundary condition is appropriate 


v o{ r — 1) = 0. (26) 

Once u -2 is determined, it is quite straightforward (although a numerical task) to deter- 
mine t>o(f)> u>o( r )i numerical results will be presented in the following section. It is worth 
noting that no, u 2 and wq, as described above, also represent an exact solution of the lin- 
earised Navier Stokes equations ((9), (10)) in the absence of the uniform flow. 

We now turn to consider the leading order terms in (17) that are dependent upon both 
r and 6 and which are, indeed, determined primarily through inviscid equations. Again, the 
linear nature of (9), (10) greatly simplifies the solution technique, this time allowing us to 
write 


iq = Ft (r) sin 9, 
Ui = «i(r)cos0, 
idi - u?i(r)sin0, 
Po = Po ( r ) cos ^ 


(27) 


6 


We may also write 


Voi ~ Coicosfl, 

U m = U 01 sin 6, (28) 

with 

= 0 -^). 

Um = -(1 + ±). (29) 

r L 


Taking terms 0(Cl Re) in the radial, azimuthal and axial momentum equations respec- 
tively, leads to 


_ dv o dV oi Vi 2uj dp 0 

V oi “j — r v o~; ( — 7 T" — ““j - ? 

dr dr r l r l dr 

dU oi 1 ^oi 1 _ 

vo—; 7^1 + v 0 = -p 0 

dr r l r r 

— dw 0 1 _ _ 

mh — \r — Ap 0? 




(30) 

(31) 

(32) 


whilst continuity leads to 

- -j-(rui) — — + AWi = 0. (33) 

7' dr r 

On this system we may impose decay of all components of the solution as r — ► oo, whilst on 
r = 1, we may have only 

tJi(r = 1) = 0. (34) 


The no-slip constraints on U] and uJ\ are therefore violated, but these may be rectified by 
the inclusion of a thin boundary layer of thickness 0(fl _ 5 Re" 5 ) on 7’ = 1. 

Specifically, we write 


and so 


giving 


V = 

(r-l)ft* Re* =0(1), 

(35) 

u, = 

ih{Y)e ie + c.c., 

(36) 

m = 

wi(Y)e ,e + c.c., 

(37) 

iu\ 

= U]YY — *Po(0> 

(38) 

ixb\ 

= W\YY + ^Fo(0) 

(39) 

U\ = - 

-Po(l)[l-e- (,+ ' |17 i, 

(40) 

uq = - 


(41) 


We see therefore, that on the 7- — 1 = (9(1) scale, intriguingly both viscous and inviscid 
effects are important simultaneously, and to a large degree, independently. In the following 
section we consider a number of numerical results arising from the results of this section, 
and go on to consider the limits A — > oo and A — > 0. 
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4 Results and large/small A behaviour 

The system (22), (25), (30)-(34) was solved using a conventional fourth order Runge-Kutta 

method. Results are shown in figure 1 for p Q (r — 1) ( ), pi(7 - = 1) (solidus), u- 2 T (v = 

1)(— • — • — ), where the dependence of these quantities with A is shown. 

We may make some further analytic progress by considering the limits of large and small 
A. Taking first A — ► oo, then by (24) 


-Y 


where 


Y = A(r - 1) = 0(1), 


i.e. 


u 2 * e , (42) 

(43) 

u 2 r| r .j -» -A. (44) 

The above implies the perturbation to the flow is confined to within a thin r — 1 = 0(A -1 ) 
boundary layer (although we must impose the restriction that A = o(fl5 Re?) in order that 
this layer remains outside the Y = 0(1) layer discussed in the previous section). After some 
algebra, it is possible to show that 


u ° — 7 T 2^ 2e Y + ^(A 3 )> 




4A 2 

,-Y 


w 0 = 

Pi “ 
Po = 


A 2 


— e 


A 3 


5- Gd" “ 2 9 ) + 0(A ' 3) ’ 

g-K)+0(A-q, 

'f 4 + ] Lf 5 -5^-5v') + 0 (A-q. 


A 


(45) 

(46) 

(47) 

(48) 


Comparison of (47) evaluated on Y = 0 with the corresponding numerical results of figure 1 
shows good agreement as A increases. Unfortunately it is not possible to compare p Q (r = 1) 
with p 0 (^ = 0) without further substantial algebra since p 0 (Y = 0) = 0(A~ 4 ). 

The alternative limit of A — > 0 is slightly more complicated, because two key radial 
lengthscales emerge. For 7’ — 0(1), we have that (24) reduces to 


1 

It 2 ~ 

r 


(49) 


which interestingly is the azimuthal velocity component corresponding to a line vortex, i.e. 
effectively the two-dimensional result. There is also a large radial scale, R = Ar = 0(1). (We 
do have the restriction that A = o(Q) in order that over the lengthscales under consideration, 
(I oo remains dominant over U 0 1 in (16)). 

For r = 0(1), the solution develops in the form 


vo - A 2 [log At)o(j’) + 0(1)] , 


(50) 
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with the leading order general solution taking the form 

vo — Ar + Br\ogr + Cr 3 -| . (51) 

r 

For R = 0(1), we must consider the first two leading terms in the solution development, 
namely 

no = A [log \v 0 (R) + v,(R) + O((log A)" 1 )] . (52) 

The complete analytic solution for v 0 {R ) and vj (R) does not appear possible, however it 
suffices to consider the limit R — ► 0, for which 

v 0 (R ) ~ A 0 R + BoRIog R + CoR 3 + (53) 

whilst vi(R) (which is forced directly by the u 2 term on the right-hand-side of (25) ) takes 
the form 


v\(R) ~ A\ R + B\ R log R + Ci 
+ ^ - ifl(log Rf. 

Boundedness constraints demand 

C=D o = D, = 0, 

whilst matching of (49) as r -» oo with (51) as R -* 0 requires 

Bo = 

4 

B = J T 

Further imposing the two boundary conditions on r = 1 yields 

- 1 1 i 1 

»o= g r- i rlogr--. 

This then leads to 


t»0 = -A log A log r + 0(A), 
pi = -1/r 2 . 


(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 


Fortunately it is possible to determine the leading-order term of pj without substantial 
algebra (and additional numerical effort), even though p\ is an order lower in log A than may 
first be expected; equation (4.19) yields the result that pj(r — 1) — » — 1 as A — ► 0, a result 
that agrees with our numerical results. However it is not possible to determine p 0 (r =1) 
without further substantial algebra and numerical effort. 

In the following section we go on to consider results for more general classes of two- 
dimensional flows involving closed streamlines, which are perturbed in some three-dimensional 
manner. 
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5 The periodicity requirement and the associated in- 
tegral conditions 

The result of Batchelor ( 1 ) may be viewed as arising from a condition of periodicity within 
the region of closed streamlines of the various physical flow quantities. It is needed because 
the Euler equations are, themselves, not sufficient to enforce periodicity; an alternative 
viewpoint is that the Euler equations are unable to capture the mean flow physics, i.e. the 
flow corresponding to zero wavenumber, a result that is little surprising. This point is equally 
important in the three-dimensional context, and is now investigated. 

We may write (9), (10) in the form 

Uo Aw + uA wo = V H -f- Re ^ V A w, (61) 

where H is the 0(e) total lead, i.e. 

H = p + u 0 • u, (62) 


and we have written the vorticity vector 

co =; coq T cw -f- 0(e 2 ). 


(63) 


Following Batchelor (1) and Blennerhassett (2) the work of this section is most efficiently 
carried out in terms of a coordinate system based on the two-dimensional , undisturbed 
streamfunction, specifically in terms of (^ 0,^,2 ). Here \j) represents the streamfunction of 
the undisturbed (u 0 ) flow, and ( is orthogonal to ?/> and z\ z remains the axial coordinate. 
If the undisturbed flow Uo is irrotational, then £ may be taken to be the standard velocity 

potential <j). The infinitesimal line element is ( — , h 2 d£, dz ) . Here h 2 is the £ coordinate 

V ?0 / 

metric, which in the case of an irrotational Uo is merely l/qo, and qo = juo|. 

Referred to the (ij>,£,z) coordinate system, the velocity vector may be written 


u = (ev cos Act, <70 + f u cos Xz, ew sin Az) + 0(e 2 ), (64) 

and the pressure as 

p = Po + f/h cos Az -f 0(e 2 ). (65) 

Referred to the (t/>,£,z) coordinate system, we may write the vorticity terms as follows 


w 0 


(0, 0, w 0 ) 


( 66 ) 


W = 


(wi sin A z, W 2 sin A z, w 3 cos A z) 

(tv^ + \uh 2 ), sin \z{ — \v — q 0 w ,/,) , 


V *2 

qo cos A 
h. 


d ,,_d_ (£ 

dxl> U 2 d( V<7o, 


(67) 
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We may then write the three components of (61) in the following form 


ULtJO 

<ll 


h.-2 


+ 


-VLU 0 


uh 


2 iff 


<7o 

Re“’ 




1 dp | 

<7o dip h 2 q'o 

1 dp i d 


1 d ( v 
h 2 d(, \<7o 
du 3 


— Xh 2 UJ 2 


+ Re 




du>3 

AUJ-\ — qo~ 


X p Re" 1 

— + -J— 
(]o h 2 


dip 

Sj(h 2 L02) - Sj l— 

dip d( \ qo 


The continuity equation, in terms of these variables is 


d . , . d ( u\ Xwh 2 

oT (M+ o7 - + = 0. 

dip dt, \<]o J qo 


( 68 ) 


(69) 

(70) 


(71) 


If we integrate each of (68)— (7 1 ) around a complete circuit in £ (lying entirely within a region 
of closed streamlines), then periodicity of the flow demands 



+ 


uh 2 

Ull’lxlr QO l/t 

Qo 


h-2 dp 
qo dip 


Re" 1 

<7o 


du 3 

~dt 


— A/l-2^2 


rfc = 0, 


(72) 



Re 1 Aw, - 7o^j | h 2 d( = 0, 
+ = °- 


h 2 dip 


(73) 

(74) 


These results are quite general, and indeed exact within the framework of the linearised 
Navier Stokes equations (9), (10). If we now direct our attention on the problem considered 
earlier in this paper, then uo is irrotational, and so in line with our comments regarding £ 
above, we can write 



Xu>2 

Qo 




= 0 , 


/ 
-Xp 


Xu)] diO 3 

dij) 


qo 


= 0, 


2 + Re- 1 4:(^)K = °- 


qo dip 

Note that the result of Batchelor (1) is retrieved from (76) by allowing A 


(75) 


(76) 

(77) 
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At this stage we can partly confirm the solution development of Section 3. Equation (77) 
above strongly suggests that oj 2 = 0 (Re), which in turn suggests that v and/or w must 
also be O(Re). This is partly confirmed by (75). Indeed, it can be shown that our solutions 
given in Section 3 do satisfy the above conditions, which clearly illustrate the importance of 
viscosity on the r — 1 = 0(1) scale. 


12 



6 Conclusions 


In this paper we have considered the effects of three-dimensionality introduced into the 
problem of a rapidly rotating circular cylinder in a uniform flow. The effects include a 
substantial (0(Re)) response in the three velocity components, in the bulk of the fluid. This 
effect is confirmed by the analysis of Section 5 in which three integral conditions were derived 
(analogous to the solitary integral condition in Batchelor’s (1) two-dimensional work). All 
the indications are that this massive response with three-dimensionality will be a generic 
feature of similar flows. It is likely that the results of this paper have repercussions for 
important practical applications, perhaps the most important of which being that of high 
lift aerofoils. 

This paper has deliberately focused on the large rotation rate problem (i.e. ft — *• oo). 
However the extremely important regime of ft = 0(1) remains to be studied. This regime 
involves a number of additional questions, perhaps the key aspect being the nature of the 
solution as the closed streamline region of the base (two-dimensional) flow is exited. It could 
well be that some form of mild discontinuity exists, for which a thin shear layer would be 
required. This whole problem woidd be a non-trivial numerical undertaking, but it is to be 
expected that the O(Re) velocity scales will persist. Equally we have implicitly assumed 
here that ft 7^ O(Re^), n positive integer, although if this condition is relaxed, it seems 
likely that the solution will be modified, but in a relatively minor manner. 

The c = 0(1) problem, at this stage, would appear to be a formidable task, quite possibly 
involving a solution of the full Navier .Stokes equations. Finally, although our study has 
concentrated on sinusoidal disturbances, in the axial direction, it is of course straightforward 
to extend our ideas to non-periodic axial disturbances using standard transform techniques. 
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Figure 1 Variation of p 0 (r — 1), P\{ r — 1) an< ^ u 2 t{ 7 ' — 1) with A. 


15 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 


pyanc reoofTiftg ouraen *q* ?hi* collection of information >s estimate to average * hour oer 'esporse, spc.taing the time *0‘ reviewing in«ruCtiO"i. sear:* r.z st aata source* 
gathering ana maintaining the data needed, and compiling ana review ng the':o!!ect‘on of rn'ormation Sena comments reaaramg this burden estimate or an, asoe:t o f this 
coneaion of information. >r.ciud>ng suggestions for reducing this ouraer. t: Washington Headouarters Services. Directorate fo' information Ope r aticns and 'es-crts, 12'S .ef’fersc- 
Davis High Aav. Suite 1204. Arlington, vfi 22202^1302. ana tc ine O'fice of Management and Budget Paperwork Reduction Proiert »C?G4-2 * 88) Washington. DC 2C5C3 


AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 

August 1993 


4. TITLE AND SUBTITLE 


3. REPORT TYPE AND DATES COVERED 

Contractor Report 


5. FUNDING NUMBERS 


THE THREE-DIMENSIONAL FLOW PAST A RAPIDLY ROTATING 
CIRCULAR CYLINDER 


6. AUTHOR(S) 

James P. Denier 
Peter W. Duck 


7, PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23681-0001 


C NAS1-19480 


WU 305-90-52-01 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


ICASE Report No. 93-52 


9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 
National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 

NASA CB-191512 
ICASE Report No. 93-52 


11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Michael F, Card 

Final Report 


12a. DISTRIBUTION /AVAILABILITY STATEMENT 

Unclassified - Unlimited 


Submitted to Quarterly J. 
Mechanics and Applied 
Mathematics 


12b. DISTRIBUTION CODE 


Subject Category 34 


13. ABSTRACT (Maximum 200 words) 

The high Reynolds number (Re) flow past a rapidly rotating circular cylinder is in- 
vestigated. The rotation rate of the cylinder is allowed to vary (slightly) along 
the axis of the cylinder, thereby provoking three-dimensional flow disturbances, 
which are shown to involve relatively massive (0(Re)) velocity perturbations to the 
flow away from the cylinder surface. 

Additionally, three integral conditions, analogous to the single condition determined 
in two dimensions by Batchelor (1), are derived, based on the condition of period- 
icity in the azimuthal direction. 


15. NUMBER OF PAGES 

19 


14. SUBJECT TERMS 

rotating cylinder, three-dimensional flow 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified 


NSN 7540-01 -280-5500 


Standard Form 298 (Rev 2-89} 

Prescribed by ANSI Std Z39-18 
29B-102 


☆ U S. GOVERNMENT PRINTING OFFICE: 1 993 - 72S-064/B6044 




















